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Abstract 

We consider stochastic equations of the prototype du(t, x) = (Au(i, x) + u(t, x) 1+f3 ) dt+ 
Ku(t, x) dWt on a smooth domain D C H d , with Dirichlet boundary condition, where 
/3, k are positive constants and {Wt, t > 0} is a one-dimensional standard Wiener 
process. We estimate the probability of finite time blowup of positive solutions, as well 
as the probability of existence of non-trivial positive global solutions. 
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1 Introduction 

Let D C R d be a bounded domain with smooth boundary dD. We consider a semilinear 
equation of the form 

du(t,x) = (Au(t,x) + G(u(t,x))) dt + Ku(t,x) dW t , t > 0, 
u(0,x) = f(x) ^0, x E D, (1) 
u(t,x) = 0, t^O, x G dD, 

where G : R — > IR+ is locally Lipschitz and satisfies 

G{z) > Cz 1+(3 for all z > 0, (2) 

C, (3 and k are given positive numbers, {W t , t > 0} is a standard one-dimensional Brownian 
motion on a stochastic basis [Tt,t ^ 0), P), and /:£)—> IR+ is of class C 2 and not 

identically zero. We assume (j2J) in sections 1 to 3 only; it is replaced by (TTI!) in section 4. 
Since we do not assume G to be Lipschitz, blowup of the solution (TjQ) in finite time can not 
be excluded, and our aim is to give estimates of the probability of blowup and conditions 
for the existence of a global solution of A (random) time T is called blowup time of u if 

lim sup \u(t, x) \ = +oo P — a.s. on {T < +oo}. 
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In the classical (deterministic) case where G(z) = z 1+/3 and k = 0, it is well-known that for 
a nonnegative / G L 2 (D), the condition 



I f\x)ip{x) dx > A 

D 



V/3 



(3) 



already implies finite time blowup of ([I]). Here Ai > is the first eigenvalue of the Laplacian 
on D, and i/j the corresponding eigenfunction normalized so that HV'IU 1 = 1. 

The existence, uniqueness and trajectorial regularity of global solutions of parabolic equa- 
tions perturbed by a time-homogenous white noise have been investigated by different meth- 
ods (see e.g. Chueshov and Vuillermot [3], Denis et al. [5], Gyongy and Rovira [9], Krylov 
[TT] . Lototski and Rozovskii [T2], Mikulevicius and Pragarauskas [33]). Several types of so- 
lutions have been proposed (see especially the last cited reference for strong solutions), and 
the regularity results show that the solution is much smoother in the space variable than for 
equations perturbed by space-dependent white noise. 

Let us recall the notions of weak and mild solutions of (1) we are going to use here. Let 
r < +oo be a stopping time. A continous ^-adapted random field u = {u(t,x), t > 0, 
x G D} is a weak solution of (1) on the interval ]0,r[ provided that, for every ip G C 2 (D) 
vanishing on dD, there holds 

/ u(t, x)f{x) dx = I f(x)<p(x)dx + / / [u(s, x)Aip(x) + G(u(s, x))(p(x)] dx ds 
J d Jd Jo Jd 



for all t G [0, r[. Let {St, t ^ 0} be the semigroup of rf-dimensional Brownian motion killed 
at the boundary of D. A continous jF r adapted random field u = {u(t, x), t ^ 0, x G D} is 
a mild solution of (1) on the interval ]0, r[ if it satisfies 



u(t,x) — Stf(x)+ / [S t - r (G(u(r, -))(x) dr + nSt- r (u(r, -))(x) dW r ] P-a.s. and x-a.e. in D 



for all t G]0,r[ (see e.g. [16], Chapter IV). We refer to [9] for background on existence of 
weak and mild solutions, and for their equivalence under local Lipschitz conditions on G. 
Let us note that the results in [9] hold for a more general class of second order differential 
operators which includes the Laplacian as a special case. The positivity of the solution of 
([1]) follows from comparison theorems (see e.g. Berge et al. |2j or Manthey and Zausinger 



Our aim in this communication is to study the blowup behaviour of u by means of a 
related random partial differential equation (see (j3J) below). In section 3 we describe the 
blowup behaviour of the solution v of this random partial differential equation in terms of 
the first eigenvalue and the first eigenfunction of the Laplace operator on D. This is done by 
solving explicitly a stochastic equation in the time variable which is obtained from the weak 
form of (|1]). The solution of this differential equation can be written in terms of integrals 




P - a.s. 
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of exponential Brownian motion with drift. The results of Dufresne [B] and Yor [T8] on 
the law of these integrals easily imply estimates for the probability of existence of a global 
solution, or of blowup in finite time of u and v. In section 4 sufficient conditions for v to be 
a global solution are given in terms of the semigroup of the Laplace operator using recent 
sharp results on its transition density. These conditions show in particular that the initial 
condition / has to be small enough in order to avoid for a given G the blowup of v, and that 
the presence of noise may help to prevent blowup. The results of section 4 can be used to 
investigate the blowup behavior of u by means of conditions (|2J) and (TTTT) . 



2 A related random partial differential equation 

In this section we investigate the random partial differential equation 

^(t,x) = Av(t,x)-—v(t,x) + e~ KWt G(e KWt v(t,x)), t > 0, x G D, 

v(0,x) = f(x), xeD, (4) 
v(t,x) = 0, xedD. 

This equation is understood trajectory wise and classical results for partial differential equa- 
tions of parabolic type apply to show existence, uniqueness and positivity of a solution up 
to eventual blowup (see e.g. Friedman [7] Chapter 7, Theorem 9). 

Proposition 1 Let u be a weak solution of (QJj. Then the function v defined by 

v{t,x) = e~ KWt u(t,x), t>0, xeD. 

solves Q). 

Remark Proposition [TJ implies in particular that (JH) possesses a strong local solution. 
Proof. Recall that Ito's formula states that {e~ KWt , t > 0} is the semimartingale given 



e -nw t = e - KW °dW s + — [ e~ KW ° ds. 

Jo 2 J 



by 

Let us write u(t, ip) = J D u(t, x)(p(x) dx. Then a weak solution of ([TJ) can be written as 

u(t, ip) = u(Q, ip) + / u(s, Aip) ds + / G(u)(s, ip) ds + k / u(s,(p)dW s . 
Jo Jo Jo 

Therefore, for ip fixed, {u(t, ip)l[ 0jT [(t) , t > 0} is again a semimartingale. By applying the 

integration by parts formula (see e.g. Klebaner [TD], Ch. 8) we get 

v(t,tp) := / v (t, x)(f(x) dx 
Jd 

k 2 



v(0,(p) + J e- KWs du(s, <p) + J u(s,<p)^-Ke- KWs dW s + —e- KWs ds 
+ [e-" w ;u{;V>)] (*), 
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where the quadratic variation is given by 

[e~ KW - } u(-, <p)] (*) = -/ K 2 e- KW °u(s } if) ds, t > 0. 

Jo 



Therefore, 



v(t,(p) = v(0,ip)+ e~ KWs (u(s,A(p) +G(u)(s,(p)) ds + k / e~ KWa u{s, if) dW t 



nt 2 pt p 

-k / e- KW °u(s, if ) dW s + ^ / e- KW/s M(s, c?s — k? 

Jo 2 y j 



v(0,(p) + 



K 



v(s, Atp) + e- KWs G(e KW v)(s, <p) - —v(s, <p) 



-kW, 



ds. 



u(s, (p) ds 



Moreover, by self-adjointness of the Laplacian, and the fact that <p(x) = for x e dD, 
v(s,Acp)= J v(s, x) A<p(x), dx = I Av{s, x)(p(x) dx = Av(s, ip). 

J D Jd 



3 An estimate of the probability of blowup 

Without loss of generality, let us assume that C = 1 in (j2J). Let ip be the eigenfunction 
corresponding to the first eigenvalue Ai of the Laplacian on D, normalized by J D ip(x) dx = 1. 
It is well-known that ip is strictly positive on D. Due to Proposition [T] we have that 



K 



v(s,Aip) - —v(s,ip) 



ds+ I e- nWs G{e KW -v){s^)ds. 
'o 



Moreover, 



and, due to (T2]) 



v(s,Aip) = -Xxv(s,ip), 



- KWs G(e KWs v(s,x))ip(x)dx > e Kf3Ws I v(s, x) 1+/3 ip(x) dx. 



D 



D 



By Jensen's inequality 

v (s, x) 1+f5 ?p(x) dx > 



D 



v(s, x)ijj(x) dx 



D 



1+/9 



and therefore 



j f v(t, j>) > - Li + y ) v(t, VO + e^v(t, 



(5) 



(6) 



(7) 



Hence v(t,ip) > I(t) for all t > 0, where /(•) solves 



d_ 
~dt 



I(t) 



Ai + y ) J(f) + e^J(t) 1 ^ , 1(0) = v(0, 



and is given by 



1(f) 



-{X 1 +K 2 /2)t 



, 0<t<T, 



with 



t := inf U > 







(8) 



It follows that I exhibits finite time blowup on the event [t < oo]. Since I ^ v(-,ip), r is an 
upper bound for the blow-up time of v(-,ip), and therefore for the blowup times of v and u. 



Remark 2 The same formula for the blow-up time of a stochastic differential equation, 
containing a stochastic integral with respect to W, has been obtained in Bandle et al pQ. 
The argument based on the first eigenvalue (and the corresponding eigenfunction) of the 
Laplace operator on D is applied there directly to u, and leads to a stochastic differential 
inequality for u(t,tp). The associated stochastic differential equation can again be solved 
explicitly by means of the Ito calculus, and the same formula as above is obtained for the 
blowup time of the solution of this equation. Both approaches are therefore equivalent, but 
the approach in p] requires a more complicated comparison theorem for stochastic differential 
inequalities. 



Let us now give an estimate for the probability of blowup in finite time of v. From ([8]) 



P[ T = +oo] = P 
= P 
= P 



ft i 

/ exp(-(Ai + k 2 /2)(3s + k(3W s ) ds < -v{0, for all t > 
Jo P 

-poo ^ 

/ exp(-(Ai + k 2 /2)(3s + k(3W s ) ds ^ -v(0, 
Jo P 



(9) 



oo i 

exp(2pW^)ds ^ -t»(0,^) _/3 



Uo P 



where Wr } ■= + W s , /j, := -(Ai + k 2 /2)/k, and (3 := k/3/2. Setting /2 = we get 

>0 



k 2 (3 2 



P[r = +oo] = P 
It follows from [18] (Chapter 6, Corollary 1.2) that 

exv(2Wp) ds 



exp 



(10) 



2Z_ 
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in distribution, where Z_-^ is a random variable with law r(— ju), i.e. P(Z_^ e ofa/) = 
T ^~^ e~ y y^~ l dy. We get therefore (see also formula 1.10.4(1) in [3]) 

,i„(0,V>)-' 3 

P[r = +00] = / h(y)dy, 
Jo 

where 

/ 2 \ 

" W yr((2A 1 + « 2 )/( K 2 /5)) 6XP I kW' 
In this way we have proved the following 

Proposition 3 The probability that the solution of (1) blows up in finite time is lower 
bounded by fi^Q^-p %) dy. 

Remark 4.1 Notice that formula 1.10.4(1) in [3] expresses the probability density function 
of J *exp(— (Ai + k 2 /2)f3s + k/3W s ) ds in terms of the Kummer functions for < 2. 

Remark 4.2 By putting k = we get v = u and, moreover, in (Q we obtain that 
P[t = +00] = or 1 according to J D f(x)ip(x) dx > or J D f(x)ip(x) dx < X 1 / 13 , which is 
a probabilistic counterpart to (j3J). 



4 Non explosion of v 

We consider again equation (Tjfl), but we assume now that k 7^ and that G : R4. — > 1R+ 
satisfies G(0) = 0, G(z)/z is increasing and 

G{z) < Az 1+I3 for all z > 0, (11) 

where A and (3 are certain positive numbers. Let {St, t > 0} again denote the semigroup of 
(i-dimensional Brownian motion killed at the boundary of D. Recall that Equation (j4]) can 
be re-written as 

v(t, x) = e~ KH ' 2 S t f{x) + f e'^-^S^r (e~ KWr G (e KM S(r, •))) (x) dr. (12) 

Jo 

We give now a sufficient condition for the existence of a global solution of (jl]). 
Theorem 5 Assume that f satisfies 

POO 

A/?/ e^We-^SrfW^dr < 1. (13) 
Jo 

Then Equation (fj]) admits a global solution v(t, x) that satisfies 

e- K ' 2t/2 Stf(x) 

< v(t, x) < r, t > 0. (14) 

(l - Af3 J^e^We-^/^Srffoodry 
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Proof. Defining 



B(t) = 1 - A/3 / e 



K/3Wr II — K^r 



we get 5(0) = 1 and 



which implies 



(t)=Ae^1e-^/ 2 ^/||^ J B 1 ^(t), 



t > 0, 



B{t) = 1 + A f e^We-^SJWiB^ir) dr. 
Jo 

Suppose now that (t,x) i— > Vt(x) is a nonnegative continuous function such that Vj(-) G 
C (£>), t > 0, and 



e - K ^ 2 S t f(x) < V t (x) < B{t)e- KH / 2 S t f{x), t > 0, x E D. 

Let t 

i?(V)(t,x) := e- K2 * /2 5 t /(x) + / e- KWr e- K2{t - r)/2 S t - r (G(e KWr V r (-))) (x) dr. 

Jo 

Then, 
R(V)(t,x) 
= e- K2t ' 2 S t f\x) + 



(15) 



e -*Wr e 



-^( t -,)/2 5 ? G(e^5(r)||e-^/ 2 g r ./l| 0O ) \ 
< e" t2 * /2 5 i /(x)+A / e ft/3 ^5 1+/3 (r)||e- K2r/2 5 r /||^e- K2( '- r)/2 ^_ r (e- K2r/2 5 r /)(x)dr 



< e-^ 2 S t f(x)+ / e-^e 



- k2 * /2 ^/(x) 



1 + A f e ^ m -B 1+p (r) ||e- K2r/2 5 r /||^ dr 
Jo 



e-^ 2 S t f{x)B{t), (16) 



where to obtain the first inequality we used the rightmost inequality in ( fT5l) and the fact 
that G(z)/z is increasing, and to obtain the second inequality we used (TTTj) . Consequently, 



Let 



e - K H/2 S j^ < R (y)(t, x) < B{t)e- K2t/2 SJ{x), t > 0, xeD. 



v t {x):=e- Kt/2 S t f{x) and < +1 (x) = R{v n ){t, x), n = 0,1,2,.... 



Letting n — ► oo yields, for t > and x £ D, 

< v(t,x) = lim u t n (x) < 5(t)e- K2 * /2 ^/(a;) < 



- 2 */ 2 SJ(x) 



1 - A(3 j*e^ w r\\e- 2 r/2S r f\\Ldr) . 



Hence, v(t,x) is a global solution of (1121) due to the monotone convergence theorem. 



Remark. If we modify (JH) and (|T2l) by replacing G(e KWr v (t, x)) by G(v (t,x)), then a global 
positive solution still exists for all / small enough, even if the inequality in (fTTj) holds only 
for z G (0, C*), where C* is some positive constant. In fact, if / satisfies 

, < C* (l- A0 J™ e- KW *\\e-* r ' 2 S r f\\£,d r y , (17) 

then Theorem still holds if we replace the factor e K/3Wr in (1131) and (|14p by the factor 
e~ nWr . We only have to verify that assuming z G (0, C*) in ffTTl) already implies the second 
inequality in (fIBl) : 



< C* (l - ^ j\-^\\e-^' 2 S r f\\^dr^ ' = for all t £ 0, 



where 

ft 



B*(t) = (l-Apj^ e-^We-^Srff^dr^J 



This yields 

B'WWe-^StfWooefrC*) (18) 

for all t > 0, since / ^ 0. 



Let us now proceed to derive a sufficient condition for (|T3l) in terms of the transition 
kernels {pt(a:, y), t > 0} of {S^, t > 0} and the first eigenvalue Ai and corresponding eigen- 
function ip. We recall the following sharp bounds for {p t (x,y), t > 0}, which we borrowed 
from Ouhabaz and Wang [T5] . 



Theorem 6 Let ip > be the first Dirichlet eigenf unction on a connected bounded C ,a - 
domain in R d ; where a > and d > 1, and /et p t (x,y) be the corresponding Dirichlet heat 
kernel. There exists a constant c > sitc/i inai, /or any t > 0, 

max < L. -/ ■'-"-><-.( sui) ^ I + r( I A A-(d+2)/2 e -(X 2 -Xi)t 



(l, It-C^)/ 2 | < e A lt gup _P^L < 1 + C (i A *)- 



where A2 > X\ are the first two Dirichlet eigenvalues. This estimate is sharp for both short 
and long times. 

The above theorem is useful in verifying condition ( 1T31) . Indeed, let the initial value / > 
be chosen so that 

f(y)<KS v i;(y), y G D, (19) 
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where rj > 1 is fixed and K > is a sufficiently small constant to be specified later on. 
Therefore S t f < KS t + v ip, and for any t > 0, 



S t f(x) < K / p t+v (x,y)ij(y)dy 



D 



\i(t+v) Pt+ri{X,y) -Xi(t 



< k(swpiP(x)) [ 

\xeD J Jl 

< K (s\xpip{x)\ I (l + c(l A (t + ■)])) 

\x£D J Jd 

= K I sup ip(x) I (e~ 



e"^""" sup - V ■ e ~^' n) ip{y)dy 



-(d+2)/2-(X 2 -\ 1 )(t+ v )\ -X t (t+r,) 



ip(y) dy 



M(t+v) + ce -Mt+v)} / ^ y )dy 
Jd 

< K(l + c)e- Xir] (supV(x) ) e~ Xlt [ ip(y)dy, 
\xtD J Jd 

which is independent of x. Since the function (t,x) i— > Stf{x) is uniformly bounded in x, 
condition ffTB"]) is satisfied provided that 



A/3 



AT(1 + c)e~ Xiri (supi{)(x)] I ip(y)dy 
\xeD J Jd 



dr e Kl3Wr-{\i+K 2 /2)/3r < l 



or 



^ r eK /3W r -(Ai+ K 2 /2)/3r 



< 



A/3 



+c) (sup^(x) ) / ip(y)dy 
\x£D J Jd 



P' 



(20) 



Notice that condition f[T7|) is satisfied if K in ([T91 is sufficiently small. 
We have proved the following 



Theorem 7 Let G satisfy / TO]) . an<i /ei D be a connected, bounded C 1,a - domain in M. d , where 
a > 0. If / TTgj) and / f^Oj) ZioZd /or some r/ > and > 0, then the solution of Equation flU| ) 
is global. 



Remarks. 1) The integral on the left of (12"U|) coincides with the corresponding integral 
in sections 2 and 3. The same type of bounds as in Section 3 can therefore be applied to 
estimate the probability of existence of a global positive solution. By means of the law of the 
iterated logarithm for W we see from (1201) that the presence of a noise may help to prevent 
blowup in finite time. 

2) If G(z) = Az 1+lS , the results of this section can be applied to the solution u of equation 
JU) since v(t,x) = e~ KWt u(t,x), t > 0, x E D. 
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